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O ! Abstract 

The LHC data on the forward heavy flavour hadron production can be a new unique 
source for estimation of intrinsic charm and bottom contributions to the proton. 
For example, we analyze the forward heavy baryon production, namely A^-baryon, 
within the soft QCD quark gluon string model and present the predictions for ob- 
servables which could be measured at the LHC. We also present some predictions 
for the D-meson production in pp collisions made within the perturbative QCD 
including the intrinsic charm in the proton that can be verified at the LHC. 



1. Introduction 



>- . There are successful phenomenological approaches to description of the soft hadron- 



nucleon, hadron-nucleus and nucleus-nucleus interactions at high energies based on 
the Regge theory and the 1/N expansion in QCD. For example, they are the quark- 
gluon string model (QGSM) [1] and the dual parton model (DPM) [2]. The main 
components of these models are the parton distributions in hadrons (PDF) and the 
^ ! fragmentation functions (FF), which describe fragmentation of quarks to hadrons. 

The PDF and FF are expressed in terms of intercepts of the Regge trajectories 
a R (0). 

^ ■ In these models the largest uncertainty in the calculations of the yields of heavy 

flavours is due to the absence of any reliable information on the transfer momentum 
t dependence of the Regge trajectories of heavy quarkonia (QQ). If teh QQ trajec- 
tories as assumed to be linear, the intercepts turn out to be low, for example, a^(0) 
is around -2 and «x(0) is around -8. As a result, the yield is very uncertain and 
decreases very rapidly with increasing quark mass. Furthermore, any knowledge of 
the t-dependence for in the region < t < M 2 q q^ and estimations of their 

intercepts become especially important for quantitative predictions. 

On the other side, to reduce the above-mentioned uncertainties, fix the Regge 
trajectories of the bottom 66-mesons and get information about quark and diquark 
fragmentation functions into heavy baryons one should compare relevant QGSM 
predictions with the LHC experimental data. Therefore the first goal of this paper 
is to present and discuss some predictions, for example, for the production of the 
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beauty A^-baryon and charmed D-meson in pp collisions. The second goal of this 
paper is to discuss the possibility of observing the so-called intrinsic quark com- 
ponents in pp collisions at the LHC energies [3]- [9] The idea of the intrinsic charm 
existence in the proton was first put forward thirty years ago by S.Brodsky with 
coauthors [10, 11]. They assumed the 5-quark state uudcc in the proton (Fig. 1). 
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Figure 1: Schematic presentation of a nucleon consisting of three valence quarks 
q v , quark- ant iquark qq and gluon sea, and pairs of the intrinsic charm (qf n qf n ) and 
intrinsic bottom quarks (q^c^). 

Later some other models were developed. One of them assumes the quasi-two-body 
state D°(uc)A+(udc) in the proton, see for example [3] and references therein. In 
[3]- [5] a probability to find the intrinsic charm in the proton was assumed to be from 
1 to 3.5 percent. The probability of the intrinsic bottom in the proton is suppressed 
by a factor ml/ml ~ 0.1 [12], where m c and rrib are the masses of the charmed 
and bottom quarks. Nevertheless, it was also shown that the intrinsic charm in the 
proton can result in a sizable contribution to the forward charmed meson production 
[13]. 

When the distribution of the intrinsic charm or bottom in the proton is similar 
to the valence quark distribution, then the production of the charmed (bottom) 
mesons or charmed (bottom) baryons in the fragmentation region is similar to the 
production of pions or nucleons. 

However, the amount of this production yield depends on the probability to find 
the intrinsic charm or bottom in the proton, but this amount looks too small. In 
this paper we continue our study of the forward heavy flavour production in pp 
collisions at LHC energies published in [14, 15] and present some estimations for the 
contribution of the intrinsic charm and bottom to the inclusive spectra of charmed 
and beauty baryons and mesons. In the next section the general formalism based 
on the QGSM for the forward hadron production in pp collisions at high energies 
is presented briefly. In the first part of the section Results and discussion, we 
present the predictions for the A c and A^ production in pp collisions at the LHC 
obtained within the QGSM without inclusion of the intrinsic charm and bottom 
in the proton. Then, in the second part of this section some estimations of the 
intrinsic charm and bottom contributions to the forward production of heavy mesons 
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and baryons in pp collisions are discussed. Finally, we give some recommendation 
for the forward LHC experiments to find the information on the quark intrinsic 
components in the proton. 



2. General formalism of the QGSM 

The general form for the invariant inclusive hadron spectrum within the QGSM 
[15, 2] is 

da _ 2E* da ~ 

= ^sd^ t =E anis)Mx ' Pt) ' {) 

where E, p are the energy and the three-momentum of the produced hadron h in the 
laboratory system (l.s.); E*, s are the energy of h and the square of the initial energy 
in the c.m.s of the pp-system; x,p t are the Feynman variable and the transverse 
momentum of h; a n is the cross section for production of the n-Pomeron chain (or 
2n quark-antiquark strings) decaying into hadrons, which are calculated within the 
quasi- "eikonal approximation" [16]; n = 1 corresponds to the left graph in Fig. 2, 
and n > 1 corresponds to the right graph in Fig. 2. 
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Figure 2: The one-cylinder graph (left) and the multicylinder graph (right) for the 
inclusive pp — > hX process. 

The function (fi n (x,p t ) is the convolution of the quark (diquark) distributions 
and the FF [15, 2, 17]: 

C(*) = F£\x + )F^\x^) + F£Xx + )F%\x_) + (2) 

2(n-l)FW(x + )4 n) (x_), 



where x± = \{\Jx 2 + x\ ± x) , and 



FW(x±) = I' dx l f^{x 1 )G T ^ h (^) . (3) 

Jx± \ X\ J 

Here r means the flavour of the valence (or sea) quark or diquark, f^ n \xi) is the 
quark distribution function depending on the longitudinal momentum fraction X\ in 
the n-Pomeron chain; G T ^h(z) = zD T ^ h (z), D T ^ h (z) is the FF of a quark (anti- 
quark) or diquark of flavour r into a hadron h (charmed or bottom hadron in our 
case), where z = x±/x\. The PDFs and FFs used by the calculations of Eqs.(2,3) 
are expressed in terms of the Regge trajectories, their intercepts «_r(0) and slopes 
ot' R (0) All the details of the calculations can be found in [15, 14, 18]. 
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3. Results and discussion 



3..1 QGSM results 

Some information on the charmonium (cc) and bottomonium (bb) Regge trajectories 
can be found from the data on the charmed and beauty baryon production in pp col- 
lisions at high energies. For example, Fig. 3 illustrates the sensitivity of the inclusive 
spectrum da/dx of the produced charmed baryons A c to different values of the Regge 
intercept aty(O). Experiment R608 [19], after measuring the decay A c — > A + 3tt, 
obtained (2.84±0.50±0.72) fib for the cross section of pp — > A° + 37r + X at \x\ > 0.5 
and -y/i = 62 GeV. The branching ratio of this decay is (2.8 ± 0.7 ± 1.1)%, therefore 
the cross section of the A c production is a(\x\ > 0.5) = (101 ± 18 ± 26) fib. The- 
oretical expectation for this cross section is 87.3 fib (cr(\x\ > 0.5) = 87.3 fib) with 
a^(0) = and 30.5 fib with a^(0) = —2.18. On the other hand, experiment R422 
[20] measured the process pp — > e~A c X at \x\ > 0.35. With a large uncertainty the 
cross section for the process was obtained to be from (26 ± 12) fib to (225 ± 9) fib 
[20]. It seems that the open circles (R608 experiment) in Fig.3 better correspond to 
our calculations (solid line) at a^(0) = 0. 

Unfortunately, there are no available data for the reaction pp — > A b X, some 
predictions for this kind of reactions are presented in [14], where it is shown that all 
the observables are very sensitive to the value of the intercept ax(0) of the T(bb) 
Regge trajectory. The upper limit of our results is reached at q;t(0) = 0, when 
this Regge trajectory is the nonlinear t-function. In fact, to measure the above- 
mentioned distributions of the process pp — > A b X one should reliably detect the A^- 
hyperon. For this purpose, we believe, the beauty baryon decays A5 — > J/^/A° — > 
fi + fi~pit~ (Fig. 4) and A b — > J/^/A° — > e + e^nn° can be used. The produced A b 
baryon can undergo the decay A& — > J/^A° with the branching ratio Br Ab ^j/^ A o = 
r Ab ->j/*Ao/r toi = (4.7 ± 2.8) • 10- 4 . The J/tf decays into fi+fi' (Br J/ ^+^_ = 
(5.93 ±0.06)%) or into e+e~ (Brj/^ e+e - = (5.94 ± 0.06)%), whereas A can decay 
into pir- (Br A o^ p7T - = (63.9 ± 05)%), or into nn° (Br A o^ n7V o = (35.8 ± 0.5)%). 

One can experimentally measure the differential cross section of this process 
da/d^pdtpdMj/y, where £ p = Ap/p is the energy loss, t p = {p in — p\) 2 is the four- 
momentum transfer, Mj/^ is the effective mass of the J/\l/-meson. In principle, the 
ATLAS and CMS detectors could detect the decay Ag ->• J/^ A fi + fi~ -n®n{it~p) 
by recording two muons and one nucleon (neutron in ATLAS - ZDC and proton in 
TOTEM ) emitted forward. However, the acceptance of the muon detection is 8° < 
&n < 172° [21]-[23], and according to our calculations the fraction of expected events 
in this region is too low. On the other hand, the electromagnetic calorimetry allows 
one to measure the dielectron pairs e + e" with acceptance about 1° < # e ( e +) < 179° 
[21]-[23]. 

In Figs. (5,6) the distributions over the proton scattering angle 9 p in the reactions 
pp — > A b X — > J/^/A X — > fi + fi~p-n~ X and over the neutron scattering angle 6 n in 
the process pp — > A b X — > J/^AqX — > e + e~ir n X are presented at the intercept 
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Figure 3: The differential cross section da/dx for the inclusive process pp — > A C X 
at a/s =62 GeV. The solid line corresponds to a^(0) = 0. The dashed curve corre- 
sponds to a^(0) = —2.18. The open circles correspond to the R608 experiment [19], 
and the dark circles correspond to the R422 experiment [20]. 

values ax(0) = (solid line) and ax(0) = —8 (dashed line). Figures 5,6 show 
a large sensitivity of these distributions to the intercepts «x = and ax = —8 
of the T(bb) Regge trajectory. One can see that the cross section da/dO Ihn for 
the nonlinear T(bb) Regge trajectory («x = 0) is larger than for the linear one 
(«x — — 8)- As is shown above for the A c production in pp collisions, the use 
of the nonlinear ip(cc) Regge trajectory results in the better correspondence of the 
calculations to the experimental data. Therefore, one can assume that all our results 
obtained at ax = are more preferable than the ones at ax = —8 and they can be 
considered as the upper limit of our QGSM calculations. The reaction pp — > A b X — )■ 
J/^A X -)■ {i + fi-(e + e-)pn-X can be measured by the TOTEM and CMS, and the 
neutron produced from the Ao decay can be be detected by the ATLAS Zero Degree 
Calorimeter (ZDC). 

In Figs. (7,8) the two-dimensional plots of the energy and scattering angle for the 
positron (electron) produced in the J/\l/ decay and the neutron produced in the A 
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Figure 4: The decay A b -)■ J/\PA° -)■ ,u + (e + )/i~(e~)p(72)7r~(7r ). 
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Figure 5: The distributions over 6 P for pp — >■ — >■ £>7r X at ^/i 
The solid (long dashed) curve corresponds to cky(O) = («t(0) = —8). 
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decay are presented. These plots correspond to ax = 0. 
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Figure 6: The distribution da/d8 n over the neutron scattering angle 9 n for pp — > 
Af,X — > e + e~n°n X at yfs = 10 TeV. The solid (long dashed) curve corresponds to 
a T (0) = (a T (0) = -8). 



One can see from Figs. (7,8) that positrons (electrons) are concentrated at the 
scattering angles 4 mrad< 8 e * + <16 mrad, whereas the neutrons are emitted mainly 
at 9 n < 0.1 mrad which could be measured with ATLAS using the ZDC (n) or with 
the TOTEM & CMS (p). The ATLAS experiment is also able to detect e+e~ by 
the electromagnetic calorimeter in the interval 1° < 6 e (e+) < 179° [21]. The ratio of 
the events presented in Fig. (7) to the total yield (let us call it the rate) is about 
4.6 percent (13.8 nb). The detailed analysis and theoretical predictions for the 
production in the forward pp collisions without the inclusion of the possible intrinsic 
bottom in the proton were presented in [14, 24]. It was shown [24] that the rate 
of the events in the reaction pp — » — > J/^> A — » e + e~ n~p at = 7 TeV 
when 500 GeV< E p < 4 TeV and the proton scattering angle 3 mrad < 9 P < 
10 mrad was about 1.33 % (40 nb). In principle, this kinematics corresponds to 
that of the TOTEM and CMS facilities at CERN. And the rate in the reaction 
pp — 7- Ab — > J/^ A — > e + e~ TT°n at the same initial energy when E n < 3 TeV 
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Figure 7: The distributions over 9 e + and E e + in the inclusive process pp — > A b X — > 
J/^A°X e + e~rm X at a/s = 4 TeV. The rate of these events is about 4.6 percent 
(13.8 nb). 

and the neutron scattering angle 6 n < 0.1 mrad is about 0.0207 (60. pb). This 
kinematics corresponds to the ATLAS and ZDC at CERN. 

Therefore, we see that the cross section of the A^-hyperon produced in pp colli- 
sions and decayed into ir~p or ir°n calculated within the QGSM without the intrinsic 
bottom in the proton can be from 60 pb to 40 nb at the different kinematics of the 
LHC experiments at y/s = 7 TeV. 

3.. 2 Intrinsic charm and beauty contribution 

Let us discuss the opportunity to find some information on the distributions of 
intrinsic charm and beauty in the proton from the analysis of the forward production 
of heavy flavour baryons in pp collisions at the LHC. First, we note that all sea quark 
distributions in the proton calculated within the QGSM give their contributions 
only to the multi-Pomeron graphs at n > 2 (Fig. 2 (right)). According to Eq. (2), 
they do not contribute to the one-Pomeron graph (Fig. 2 (left)). The one-Pomeron 
graph results in the main contribution at large values of x because a n decreases very 
fast when n increases [16] . Therefore the sea charm and beauty quark distributions 
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Figure 8: The neutron energy distribution as a function of 9 n for the process pp — > 
A b X -> e + e-7T°n X at y/a = 10 TeV. 

result in very small contributions to the inclusive spectra of the charmed and beauty 
hadrons at not large p t and different values of x [25, 15, 14]. For the forward heavy 
flavour hadron production this contribution is smaller and we neglect it. 

These sea charm and beauty quark distributions greatly differ from the distri- 
butions of the intrinsic charm and bottom quarks in the proton, which, according 
to the assumption of [10, 3], should behave like the valence quark distributions. 
Therefore, if one wants to estimate the contribution of the intrinsic charm and bot- 
tom quarks within the QGSM scheme, one should include their distributions in the 
calculation of the one-Pomeron graph (n — 1). The procedure can increase the rates 
of events for both reactions pp — > A b X — > e + e~p-K~X and pp — > A b X — > e + e~nn°X 
when the final proton or neutron is emitted in the forward direction because the 
one-Pomeron graph (Fig. 2 (left)) makes the major contribution to the spectrum 
in this kinematics. 

Assuming the existence of the intrinsic bb-pa.ii in the proton, as a pair of the 
valence quark-antiquark , with some non-zero probability w bb , one can estimate en- 
hancement in the forward A;, pp-production at the LHC The expected enhancement 
will be about the ratio of the differential cross sections %z(pp ->nl) l%{PP ~^ AbX) 
at very large x multiplied by w bb , Although it is suppressed in comparison with the 
intrinsic charm probability w cS by a factor w? c jm\ ~ 0.1 [12], it can not be neglected 
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nevertheless. 

Calculating these spectra within the QGSM [17, 14, 15] at large x and assuming 
w b i ~ 0.3% [5, 12] one can get that the yield of A b produced hadronically in the 
forward direction can increase by a factor 3-5 times due to the intrinsic bottom 
quark contribution. It means that the cross section of the forward production of 
Ab in PP collisions at LHC energies, which decays into e + e~ir~p or e + e~7r°n, can 
reach a few hundred nb for TOTEM and CMS and few hundred pb for ATLAS. Our 
estimations show that the yield of the forward charmed A+-hyperon production can 
be increased by a factor of 10 due to the intrinsic charm quarks. Therefore, the 
reaction pp — > A^X — > A°7r + A — > n7T°TT + X can also be measured at the LHC when 
the neutron is emitted in the forward direction. The neutron can be measured by 
the ZDC and 7r + -meson can be detected by the hadron calorimeter. We presented 




Figure 9: The D + D distributions over the pseudo-rapidity rj in pp — > (D + D )X 
at v/i = 7 TeV and 10 < p t < 25 GeV/c. 

the qualitative estimations for the contributions of the intrinsic beauty and charm 
to the forward Ab and A c production at LHC. These spectra were calculated within 
the nonperturbative QGSM in which the PDF do not include the intrinsic charm or 
beauty contributions. It is not so easy to take into account the intrinsic charm (IC) 
contribution at the PDF used in the QGSM [1, 17]. 

However, there are the PDF used in the perturbative QCD calculations which 
include the IC contribution in the proton [3]- [5]. The probability distribution for the 
5-quark state (uudcc) in the light-cone description of the proton was first calculated 
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Figure 10: The double D° distributions over the pseudo-rapidity rj in pp — > (Dq + 
D )X and at p t > 25 GeV/c. 



in [10]. It has the following form [4] 

^ = f c ( x ) = f-( x ) = J\f x 2 x ((1 -x)(l + Wx + x 2 ) + 6x(l + x)ln(x)\ , (4) 
ax 1 ' 

where M is the normalization constant. One can see from Eq.(4) that the IC dis- 
tribution has some enhancement at large x and vanishes at x — 1. As is shown 
in [5] this enhancement starts at x > 0.2 that can result in similar enhancement 
in the inclusive spectra of the open charm at large rapidities y or pseudorapidities 
rj and transverse momenta pt- We calculated the IC contribution to the inclusive 
spectra of the D-mesons produced in pp collisions at y/s =7 TeV within the pertur- 
bative QCD. In Fig. 9 the inclusive spectrum of single Z)°-mesons is presented as a 
function of the pseudorapidity rj at 10 GeV/c< p t <25 GeV/c. Calculating these 
spectra within PYTHIA8 we used the PDF both for the CTEQ66 without the IC 
(the dashed blue distributions in Fig. 9) and the CTEQ66c including the IC with the 
probability about 3.5 % (the solid red distributions in Fig. (9)) at Q 2 = ml =1.69 
GeV 2 [5]. Figure 9 shows that some enhancement due to the IC can be visible at 
large rj. Its amount increases when p t grows and, for example, the inclusion of the 
IC increases the spectrum by a factor of 2 at rj =4.5. Similar effect was predicted 
in [26]. In Fig. 10 the inclusive spectrum of the D°-meson produced in the process 
pp — > D°D°X at the rapidity interval 3< y <4.5 and -y/i =7 TeV is presented as a 
function of pt, when two D° are produced. The enhancement of the spectrum (the 
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excess of the solid histogram in comparison to the dashed one) at 7< p t <10 GeV/c 
is not more than 30 %. The predictions presented in Figs. (9,10) can be verified by 
the LHCb experiment at CERN because their facility is able to measure inclusive 
spectra of D-mesons at rj < 4.5. 

4. Conclusion 

We analyzed production of charmed and beauty baryons in proton-proton collisions 
at high energies within the soft QCD quark-gluon string model. This approach 
described rather satisfactorily the charmed baryon production in pp collisions [15, 14] 
and allowed us to apply the QGSM to beauty baryon production. We focus mainly 
on the analysis of the forward A& production in pp collisions at LHC energies. We 
present the predictions for the reaction pp — > A^X — > e + e~mr X, which can be 
studied in the ATLAS experiment using the ZDC, and for the process pp — > A^X — > 
e + e~pit~ X , which could be reliably studed at the TOTEM and CMS within their 
common proposal for Diffractive and Forward Physics at the LHC [22, 23]. We 
show that beauty A fe -hyperons can in principle be detected in LHC experiments via 
registration of their decay products A& — > e + e~p7i~ and A5 — > e + e~rm Q . 

We would like to stress that any data on forward production in the process pp — > 
AbX could be important for determination of ax(0) of the T(bb) Regge trajectory. 
Our predictions with the nonlinear Regge trajectories of cc and bb mesons can be 
considered as the upper limit of the QGSM calculations without the intrinsic charm 
and bottom in the proton. 

The inclusion of the intrinsic bottom or/and charm in the proton can increase the 
yield of the relevant heavy flavour baryons by a factor of 3 to 10. In particular, we 
considered a possibility of measuring the reaction pp — > A+X — > A°n + X — > n-K ir + X 
using the ATLAS (ZDC). This measurement can provide information on the intrinsic 
charm in the proton, the probability of which is estimated to be a factor of 10 higher 
than the one for the intrinsic bottom in the proton. Finally, it is worth noticing that 
any reliable non-observation of this enhancement in the experiments at the LHC can 
severely constrain the intrinsic heavy quark hypothesis. 

Our calculations of the charmed meson production in pp collisions within the 
MC generator PYTHIA8 and the PDF including the intrinsic charm showed the 
following. We found that the contribution of the intrinsic charm in the proton 
could be studied in the production of D-mesons in pp collisions at the LHC. The IC 
contribution for the single D°-meson production can be sizable, it is about 100 % 
at large rapidities 3 < y < 4.5 and large transverse momenta 10 < p t < 25 GeV/c. 
For the double D° production this contribution is not larger than 30 % at pt > 5 
GeV/c and 3 < y < 4.5. These IC contributions for the single and double D-meson 
production were obtained with the probability of the intrinsic charm taking to be 
w c5 =3.5 % [5], and they will decrease by a factor of 3 when w c5 ~ 1 %. Therefore, 
this value can be verified experimentally at LHCb. 
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The presented predictions could stimulate measurement of the single and double 
D-meson production in pp collisions at the CERN LHCb experiment in the kinematic 
region mentioned above to observe a possible signal for the intrinsic charm. The 
intrinsic beauty in the proton is suppressed by a factor of 10, therefore its signal in 
the inclusive spectra of 5-mesons will probably be very weak. 
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